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Abstract 

Using ab-initio electronic structure methods we have investigated the ground state structures 
and properties of neutral and charged SiHm(ni=l-4) and Si2H„(n=l-6) clusters which are produced 
in the plasma enhanced chemical vapor deposition process used in the preparation of hydrogenated 
amorphous silicon(a-Si:H). Our results show that charging a neutral cluster distorts it and the 
distortion mainly occurs through the orientation of Si-H bond. We attribute structural changes 
in the charged clusters to electrostatic repulsion between the bonded and non-bonded electrons. 
We find that in addition to the usual Si-H bond, hydrogen atoms form Si-H-H and Si-H-Si bonds 
in some clusters. The vibrations of Si-H, Si-Si, Si-H-Si bond stretching modes show that the 
frequencies are shifted significantly upon charging. The frequency shifts in the charged clusters are 
consistent with their bond length variations. We discuss the fragmentation pathways of silane into 
binary products and the role of fragmented silane radicals in the cluster formation and a-Si:H film 
deposition process. 

PACS numbers: 73.22.-f, 31.15.Ar, 36.40.Wa, 36.40.Mr 
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I. INTRODUCTION 



Silane(SiH4) is used in the plasma enhanced chemical vapor deposition(PECVD) process 
to produce hydrogenated amorphous silicon (a-Si:H) which has potential applications in solar 
cells and flat panel devices. Before a-Si:H film deposition takes place, the gaseous silane 
molecules are decomposed either by using electric field [l| or by thermal energy The 
decomposed species live for a short time in plasma environment and contribute to a-Si:H film 
growth. In the plasma, the decomposed species exist also in charged states. Experimental 



studies 



. in tnepiasma, tne 



Id) llll have shown the presence of various neutral and charged 



species in the plasma such as Si, SIHq and SioHgand SiHs^, Si2H5+ and Sia;H+ (x=4-10). 



Also the experimental investigations! 



show that the decomposed species 



have a great influence on the properties of a-Si:H films. Therefore, investigating the structure 
and properties of neutral and charged silane and disilane derivatives can help in gaining a 
better understanding of a-Si:H film growth in the PECVD process. With this motivation 
we have carried out ab-initio electronic structure calculations to study (1) the ground state 
structure, (2) vibrational frequencies, (3) fragmentation and (4) clustering process of neutral 
and charged SiHm(m=l-4) and Si2H„(n=l-6) clusters. 

In the recent past, several electronic structure calculations have been presented on neutral 

19, 3- 
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and charged silane and disilane derivatives at various levels of accuracv|15l. 
The ground and excited states of SiH have been investigated using multireference con- 
figuration interaction approach jl5|. Using configuration interaction and coupled cluster 
method with single and double excitations, a detailed vibrational analysis on SiHs" has been 
performed 3]. Coupled cluster calculations have been carried out to study the ground state 
structures and vibrational properties of SiH„ (n = 2-4) clusters [t^. Grev and Schaefer 18 1 
applied coupled cluster method to identify the local minima of Si2H2. Gupte and Prasad IS | 
reported ground state structures and vibrational frequencies of Si„Hm (n=l,2 and m = 1-6) 
clusters using non-orthogonal tight binding molecular dynamics. Recently, the potential 
energy surface of Si2H~ has been explored |2fll| using coupled cluster and configuration inter- 
action methods. However, ab-initio calculations on neutral and charged silane and disilane 
derivatives using the same methodology are still lacking. We have performed density func- 
tional theory calculations to understand the charge induced effects on the structure and 
properties of silane and disilane derivatives. 
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The fragmentation and the clustering reactions of silane derivatives in the PECVD process 
are not yet well understood. Phenomenological model calculations have been carried out to 
understand the formation of various clusters during silane discharge and their role on film 
deposition processj^, 0]. Using semi-empirical method the collisions between Si^ and H2 
over a range of collision energies have been studiedj^l- Sosa and Lee|23j performed density 
functional theory calculations to study the insertion reactions of SiH2, SiHF and SiF2 into 



H2 and found that fluorine substitution increases t 
theoretical investigations, Raghavachari 
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le reaction barrier heights. In a series of 



27( 1 showed that the sequential clustering 



reactions of silane with Si^, SiH"^, SiH2''' and SiHa^ end with different products such as 
Si2H4"'", Si2H5+, SisHs^, SisHs^, Si4H7+, SisHio^ and SisHn"*". We have also investigated 
the fragmentation of silane and the role of the fragmented radicals in the beginning steps of 
the clustering process. 

Our ab-initio calculations on neutral and charged SiHm(m=l-4) and Si2H„(n=l-6) clus- 
ters show several interesting findings about their ground state structures and properties. We 
find that in addition to the usual Si-H bond, hydrogen atoms in some clusters form Si-H-H 
and Si-H-Si bonds. Charging the cluster affects the Si-H bond orientations and Si-H, Si-H-Si 
and Si-Si bond lengths. We explain the bond orientations and bond length variations of the 
charged clusters using molecular orbital, charge density and population analysis 

We have done vibrational analysis of the neutral and charged clusters and find that the 
frequencies of stretching modes of Si-H, Si-Si and H-Si-H bonds shift considerably upon 
charging. These frequency shifts are consistent with the bond length variations of the 
charged clusters. The fragmentation studies on neutral and charged silane show that the 
major radicals formed are SiH2, SiH2^, SiHa"*" and SiHs^. These radicals undergo further 
clustering processes which lead to disilane derivatives. We also discuss the possible role of 
these radicals in a-Si:H film deposition process. 

The plan of the paper is as follows. In section II we give computational details of the 
present work. In section III we present and discuss in detail the ground state geometries of 
neutral and charged silane and disilane derivatives. Vibrational analysis of these clusters is 
discussed in section IV. In section V we discuss the fragmentation of silane and the formation 
of bigger hydrogenated silicon clusters. Finally, in section VI we give our conclusions. 
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II. COMPUTATIONAL DETAILS 

We have used the Car-Parrinello molecular dynamics(CPMD) Q^Sl and Gaussian98 
package for our study. Using the CPMD method Isil-ls^ with simulated annealing optimiza- 
tion we have obtained the ground state structures of the neutral clusters. The pseudopoten- 



tials for si 
technique 



icon and hydrogen have been generated using the Bachelet, Hamann and Schliiter 



33l |. The local density approximation(LDA) of the density functional theory has 



been used with the Ceper 
by Per dew and Zunger 



ey-Alder 



3J] exchange-correlation energy functional parameterized 



The wave functions were expanded in a plane wave basis with 12 
Rydberg energy cut-off and k=0 point was used for Brillouin zone sampling. During simu- 
lation volume of the system was kept constant and to avoid interaction between the clusters 
a big fee supercell with side length of 35 a.u. was used. To perform simulated annealing, 
the system was taken to high temperatures (1200K in the steps of 300K), equilibrated for a 
long time (about 16,000 steps) and then slowly cooled down (in the steps of 50K) to 300K. 
Below this temperature the steepest descent optimization was found to be more efficient 
to obtain the ground state geometry. The desired temperature was achieved by rescaling 
atomic velocities and the atoms were moved according to the velocity Verlet algorithm [36| 
with a time step of 5 a.u. The fictitious mass of the electron was taken to be 200 a.u. 
Our SiH„(n = 2-4) and Si2H2 ground state structures obtained using the CPMD method 
are in good agreement with those found earlier Q| using ab-initio quantum chemical 
calculations. 

For obtaining the ground state structures of the charged clustersjs^] we have used Gaus- 
sian98 package. We have taken the structures of neutral clusters obtained from the CPMD 
calculations and re-optimized the neutral and charged clusters using Gaussian98 package by 
Berny's optimization technique . The geometry optimization has been done with Becke's 
3 parameter hybrid functional(B3LYP) for exchange- correlation functional and employing 
6-31 Ig** basis set HQ- The vibrational frequencies of the neutral and charged clusters are 
calculated with the same functional and basis set. We see from Table I that the ionization 
potentials (IP) of the clusters obtained from Gaussian calculations are in good agreement 



with the experimental values 



40 
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l42j . The IP values of Si and Si2 have been calculated 
with various multiplicities to take account of the orbital degeneracy. We find that the ground 
state of Si, Si"*", Si2 and Si2 have multiplicities 3, 2, 3 and 4 respectively. The IP values of 
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Si and Si2 reported in Table I have been obtained using the ground state multiphcities. Our 
calculated umbrella mode frequency 894 cm~^ of SiHa^ is in excellent agreement with the 
experimental value of 880 cm~^. 



III. GROUND STATE STRUCTURES 

A. Charge induced structural modifications 

In this section we discuss the charge induced effects on the structures of SiHm(m=l-4) and 
Si2Hn(n=l-6) clusters. In Fig. 1 and 2 we have shown the ground state structures of neutral 
and charged silane and disilane derivatives obtained using Gaussian98 package employing 
B3LYP/6-311G**. From Fig. 1 we see that the structure of SiH2 is asymmetrical top whose 
bond angle opens up when positively charged and does not change much when negatively 
charged. The structure of Sills is pyramidal which becomes planar when positively charged 
and Si-H bonds further bend when negatively charged. 

The ground state structure of SiH4 has tetrahedral symmetry. This symmetry is destroyed 



upon charging j41| 



44L |45|and we see in Fig. 1 that SiH4+ and SiH4 have very different 



structures. This symmetry breaking has been attributed to the Jahn- Teller effect |45 
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48|.The structure of SiH4"'" obtained by our calculation is in agreement with the structure 



amed py our caicuu 
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obtained by earlier calculations |42l, |^ |^ Our ground state of SiHs^ is in good 

agreement with the result obtained by Boo and Lee0|. The structure of SiHs"*" is more like 
a complex of SiHs^ and II2 and that of SiHs" is more like a complex of SiH4~ and H. 

Interestingly, two H atoms in SiH4"'' and SiHs^ are over-coordinated and form Si-H-H 
bond. This Si-H-H bond is formed between three atoms by two electrons and is an example 



of multi-center bond 55[ . A similar type of three center- two electron bond, C-H-H bond, has 
been observed 0. Isfi] in CH5+. It was predicted that due to the three center - two electron 



54 



55, that is, these molecules do not 



bond, SiHs"*" and CHs"*" have fiuxional behavior 
have definite structure and bonding arrangement which makes it difficult to characterize 
them using spectroscopic techniques. These results show that charging clusters not only 
changes the structure but also can give rise to the formation of unusual bonds. 

The over-coordination of H atom has also been seen in neutral clusters such as Si2H, 
Si2H2 and Si2H3. We see from Fig. 2 that hydrogen atom is over-coordinated in these 
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clusters by connecting both silicon atoms and forms a Si-H-Si bridge bond. Note that 
Si-H-Si bridge bond is a multi-centered bond and is different from Si-H bond, which is a 
two-centered bondjs^. The Si-H-Si bridge bonded frame of the clusters does not change 
when the clusters are charged. In the case of negatively charged Si2H2, hydrogen atoms are 
slightly pushed away from each other without changing each Si-H-Si frame. 

Now we discuss the charge induced effects on the structure of disilane derivatives. From 
Fig. 2 we see that Si-H bonds of Si2H3 are tilted away from Si-Si bond when positively 
charged and toward Si-Si bond when negatively charged. Si2H4 is nearly planar and results 
in a perfectly planar structure when positively charged. On the other hand, Si-H bonds of 
Si2H4 tilt significantly out of plane when negatively charged. Si2H5 consist of SiHs and SiH2 
units connected by Si-Si bond. We see from the figure that H atoms of SiH2 unit in Si2H5+ 
move slightly upward while those in Si2H5~ are further pushed down. From the figure we 
see that Si2H6 consist of two SiHs units connected by Si-Si bond. Upon charging there is 
no significant change in the overall structure, but each SiHs units of Si2H6'*' and Si2H6~ 
distort as those of SiHs^ and SiHs^ clusters. The ground states of neutral and charged 
silane and disilane derivatives show that the structural modifications occur mainly through 
the bending of Si-H bonds. 

The charge induced structural distortions of silane and disilane derivatives can be under- 
stood from the nature of valance electron density distribution. First, we discuss the nature 
of valance electron distribution and its influence on charged SiH2 structure. Then, using 
the same arguments we explain the structural modifications of other charged clusters. The 
valance electron density distribution of neutral and charged SiH2 are shown in Fig. 3. From 
the figure we see that the electron density of SiH2 is highly localized near H atoms and also 
on top of silicon atom. We call the electron density between Si and H as bonded while on top 
of Si as non-bonded. This non-bonded electron distribution is in some sense, similar to the 
Lewis dot 57| of unpaired electrons. The density of bonded electrons is highly localized near 
H atoms due to the fact that hydrogen is more electronegative j3] and withdraws electron 
from silicon atom. We see from the figure that the non-bonded electron density of SiH2, 
on top Si atom, decreases when positively charged and increases when negatively charged. 
Thus in Fig. 1 for SiH2 the non-bonded electron density is schematically represented by a 
half filled circle on Si. For SiH2"^ it is shown by an open circle indicating the decrease in 
the non-bonded electron density while for SiH2~, it is shown by a filled circle indicating the 
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increase in the non-bonded electron density. 

Now we discuss how the structural modifications of charged SiH2 can be understood in 
terms of the bonded and non-bonded electron density. Since the ground state of SiH2 is 
not a linear geometry, the non-bonded electron density would exert a downward force on 
the bonded electron density. Since the non-bonded electron density of SiH2'*' is lower than 
that of neutral SiH2, its downward push on H atoms of SiH2''' is also lower than that of 
neutral SiH2 which results in the opening of H-Si-H bond angle. From this argument one 
would expect that SiH2'^"'' should have further opening in the bond angle. Indeed, we find 
that the optimization on SiH2'''^ results in an almost linear geometry with a larger bond 
angle. In contrast, in SiH2~ there is increment in the non-bonded electron density on top 
of silicon atom. Thus one would expect that the bond angle should get smaller than that of 
the neutral cluster due to the excess push on H atoms. But the reduction of H-Si-H angle 
is opposed by the increased H-H repulsion. The H-H repulsion in SiH2~ increases because 
charges on H increase when SiH2 is negatively charged |59l|. 

The structural modifications of other charged clusters can be explained using the same 
arguments given in the previous paragraph. To understand the nature of valance electron 
distribution, we have done charge density analysis of all these clusters. We find that the 
non-bonded electron density distribution decreases when positively charged and increases 
when negatively charged. From the charge density analysis we have also located the non- 
bonded electron density in these clusters. As in the case of neutral and charged SiH2 
clusters, the non-bonded electron density is schematically represented by half, open and filled 
circles for the neutral, positive and negative clusters in Fig. 1 and 2. From this schematic 
representation it is easy to understand the structural modifications in the charged clusters. 
Hydrogen atoms bonded to Si are attracted to the nearby open circle while they are repelled 
away from the nearby filled circle. The attraction of H to an open circle and the repulsion of 
H to a filled circle are due to the electrostatic interaction as explained in the case of charged 
SiH2 clusters. 

B. Charge induced bond length variations 

We find that charging the clusters not only changes Si-H bond orientations but also Si-H, 
Si-H-Si and Si-Si bond lengths. In Fig. 4(a), 4(b) and 4(c), we have shown Si-H, Si-H-Si 
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and Si-Si bond lengths of the neutral and charged clusters as a function of cluster size. We 
see from the figure that Si-H bond length of a positively charged cluster is shorter than that 
of the corresponding neutral cluster while it is longer for the negatively charged cluster. On 
the other hand, Si-H-Si and Si-Si bond lengths of both positive as well as negative charged 
cluster are longer than that of the neutral cluster. However, Si2H3+ is an exception where 
Si-Si and Si-H-Si bond lengths are shorter than those of the neutral cluster. /^From Fig. 
4(c) we see that Si-Si bond length of a neutral Si2H„ cluster increases as a function of H 
and reaches maximum for Si2H6, but the bond length of Si2H4 is shorter than that of its 
neighboring Si2H3 and Si2H5 clusters. 

To understand the Si-Si and Si-H bond length variations of charged silane and disilane 
derivatives we first consider the bond length variations of Si2 and SiH dimers. Si2 has bond 
length of 2.144A, the bond expands in both positive as well as negative charged state to 
2.308 and 2.453A respectively. On the other hand, depending on the charged state the 
bond length of SiH dimer decreases or increases. SiH dimer has bond length of 1.536A, 
which shortens to 1.518A when positively charged and expands to 1.566A when negatively 
charged. The different trends seen in the bond lengths of charged Si2 and SiH dimers can 
be understood with the results obtained from population analysis. The population analysis 
on charged dimers are done by keeping the bond distance same as those of neutral dimers. 
The population analysis on Si2 shows that each atom of the neutral dimer has zero charge, 
whereas each atom of Si2'*' has +0.5e and that of Si2~ has -0.5e. The atoms of charged Si2 
are equally charged and repel each other which results in bond length expansion in both 
positive and negative charged state. The population analysis on neutral SiH dimer shows 
that a net charge of +0.35e on Si and -0.35e on H atom. The charge transfer from silicon 
to hydrogen is the indication of polar covalentjool bond which is consistent considering the 
fact that electronegativity of silicon is smaller than hydrogen [ssj. The population analysis 
on SiH+ shows that the silicon atom in the dimer holds a charge of +1.36e and hydrogen 
-0.36e. This means that the lost electron is primarily from Si atom in the SiH dimer. Hence 
the positively charged silicon pulls the H atom closer and reduces the bond length. In the 
negative charged state. Si atom carries a charge of -0.68e and H atom, a charge of -0.32e 
and this shows that the extra electron is accumulated primarily on silicon atom. This extra 
electron density on Si of SiH~ repels the H atom and results in the bond length expansion. 
Our population analysis on neutral and charged SiH dimer shows that Si atom holds the 
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excess charge in SiH"*" and SiH~. This can be understood by considering the energies required 
to add or remove an electron from Si and H atoms. Si atom has smaller ionization potential 
( IP ) than H atom while it has larger electron affinity ( EA ) than H atom 61[. Therefore 
adding or removing an electron from Si is easier with respect to H atom. 

The general trends of Si-Si and Si-H bond length variations of charged clusters are con- 
sistent with the bond length variations of Si2 and SiH dimers when charged. As mentioned 
earlier, Si-H-Si bond expands under both the charged state. In Si-H-Si unit, two silicon 
atoms are also bonded and they expand irrespective of the charged state. As a result of 
Si-Si expansion, Si-H-Si bridge bond expands under positive as well as negative charged 
state. 

The shorter Si-Si bond observed in Si2H4 is due to the formation of a tt bond between 
silicon atoms. The tt bond has a shorter bond length when compared to a cr bondjo^. 
The highest occupied molecular orbital(HOMO) of Si2H4, Si2H3 and Si2H3"'" are shown in 
Fig. 5. From the figure we see that the HOMO of Si2H4 is spread above and below the 
line connecting the two Si atoms. This indicates that the HOMO of Si2H4 is formed due 
to the side overlap between valance p orbitals of silicon atoms which is a characteristic of 

n 

TT bonding[6^. Si2H3+ has shorter Si-Si and Si-H-Si bonds compared to the corresponding 
neutral cluster. From Fig. 5 we see that the HOMO of Si2H3+ is slightly spread above and 
below the Si-Si bond while such a spread is absent in neutral Si2H3. This means that HOMO 
of Si2H3^ has gained slightly tt character compared to that of neutral Si2H3 and hence Si-Si 
bond is shorter than that of the neutral cluster. Due to the reduction in Si-Si bond length, 
Si-H-Si bond length of Si2H3+ is also reduced from that of the neutral cluster. 



IV. VIBRATIONAL ANALYSIS 



In this section we discuss the charge induced effects on vibrational frequencies of the 
clusters. Also we identify the modes corresponding to internal rotation of subunits in the 
clusters. The vibrational frequencies and force constants of silane and disilane derivatives 
are given in Tables II and HI. Our calculated vibrational frequencies 922, 980, 2235 and 
2244cm~^ for SiH4 are in good ag .ee,„e,.w.K .hose of e^pert-ne,... Obtained va.„esH 
of 914, 953, 2189 and 2267cm-^ 

We notice from the tables that SiH4+, SiHs"*", Si2H3+, Si2H4~, neutral and charged clusters 
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of Si2H5 and Si2H6 have one mode with very low frequency and nearly zero force constant. 
A small frequency and force constant for a mode imply that during the vibrational motion 
the atoms of the cluster move in a nearly constant potential energy surface. This ma y b e an 
indication that the mode corresponds to either transition state or internal rotationj^, Qi 



65|,|66l of the cluster. We looked into the vibrational motion corresponding to low frequency 
modes of the clusters. We find that the low frequency modes are associated with the internal 
rotation of subunits about some symmetry axis of the clusters. 

One vibrational frequency of SiHs corresponds to the umbrella mode. The umbrella mode 
frequency of SiHs is 760cm~^ which increases for both SiHs^ and SiHs^ to 839 and 894 cm~^ 
respectively. The inversion barriers for SiHa and SiHa^ are 0.22 and 1.24 eV respectively. 
Because of the increased barrier height, the frequency of SiHs" is larger than that of SiB.^. 
Although SiHs"*" is planar without inversion barrier, the frequency and force constant of the 
umbrella mode are larger than those of the neutral and negative clusters. This is due to the 
fact that the potential energy surface of SiHa^ about the center of inversion has different 
shape from those of SiHa and SiHs^. At the inversion center, SiHs and SiHa" have barriers 
while SiHa"'" has potential energy minimum as shown in Fig. 6. 

We find that charging the cluster shifts the stretching frequency of Si-H, Si-Si and H-Si-H 
bonds. From Table II and III we see that the trend of frequency variations in Si-H, Si-Si 
and Si-H-Si stretch modes of charged clusters is opposite to the trend of their bond length 
variations which is shown in Fig. 4. For example, Si-H stretch frequency of SiH dimer is 
2013 cm~^ which increases to 2137 cm^^ for SiH+ and decreases to 1827 cm~^ for SiH~. 
From Fig. 4 we see that the bond length is shorter for SiH"*" while it is longer for SiH~ 
from that of the neutral SiH dimer. Such a relationship between Si-H bond length variation 
and frequency variation seen in SiH dimer holds true for Si-H, Si-Si and Si-H-Si bonds of 
all the clusters. In general, longer bond length reduces the frequency and vice versa. This 
relationship is reasonable since expansion of the bond means that the bond gets weaker and 
hence the corresponding stretch mode has lower frequency. On the other hand shortening of 
the bond length means that the bond gets stronger and hence the stretch mode has larger 
frequency. 
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V. FRAGMENTATION OF SILANE AND THE FORMATION OF BIGGER 
CLUSTERS 



A. Fragmentation of neutral and charged silane 

For the fragmentation studies, we have considered various possible fragmentation path- 
ways of silane into binary products. Also the charged clusters are included in our investi- 
gation as they exist in PECVD plasma. The fragmentation energy of a cluster has been 
calculated by taking the difference between total energy of reactant with those of product 
clusters. Among the possible fragmentation channels of neutral and charged silane the low 
energy channels are 

SiHi S1H2 + i/2(+2.57e1/) (1) 
SiH^ SiH+ + if2(+0.60eK) (2) 
SiHl SiH^ + H{+1.33eV) (3) 
where the number in the parenthesis is the fragmentation energy. 

As we see from the above reactions, production of hydrogen in molecular form is favored in 
SiH4 and SiH4"'" while atomic hydrogen is favored in SiH4~. Our low energy fragmentation 
channel of SiH4 agrees with the experimental studyjsj. There is another process where 
hydrogen atom is attached with SiH4+ and forms SiHs^, which undergoes fragmentation as 
given below. 

SiH+ + H ^ SiH+ SiH+ + H2{+0MeV) (4) 

This process is quite possible since a large amount of hydrogen is available in the PECVD 

plasma and we see from the fragmentation energies of SiH4"'" and SiHs"*" that the latter is 

more stable to the fragmentation of H2. We find that the energy required to break H2 

from SiHs" is negative and is about -0.24 eV. This indicates that SiHs" is highly unstable. 

Therefore, SiHs" is not considered in our fragmentation study. 

It has been observed that deposition of SiH^ on the film improves the material quality 

n 

whereas deposition of SiH and SiH2 degrades the material quality 1]. Therefore, it is impor- 
tant to understand the production of SiHs and its role on the film deposition process. We see 
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from equations (3) and (4) that SiHa" and SiHa"*" are produced in the fragmentation process 
after attaching an electron with SiH4 and hydrogen atom with SiH4"'". This indicates that 
the production of SiHa can be improved by controlhng the electron and H content in PECVD 
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plasma. It is interesting to compare these processes with the experimental investigation 
on argon/hydrogen expanding thermal plasma. The investigation suggests that a significant 
amount of SiHs is produced by electron collision and hydrogen collision process. 

The structure and dipole moment of SiHs may play a role on improving a:Si-H film 
quality. As we have mentioned, SiHa has low inversion barriers of 0.22 eV. This implies 
that SiHs can transform to planar structure from its ground state pyramidal structure by 
crossing the barrier. The barrier crossing can easily occur in the hot plasma environment. 
When the cluster transforms into the planar structure, its dipole moment changes by a 
considerable amount. The dipole moment of pyramidal SiHs is 0.16D and that of planar SiHa 
is zero. Because of its zero dipole moment the planar SiHa does not get easily physisorbed 
on the film. This low physisorbtion behavior of planar SiHq cluster is similar to that of 

n 

CH3 in hydrogenated amorphous carbon film deposition |68l|. Since the planar SiHs has low 
physisorbtion property, it can migrate on the surface until it gets absorbed by a highly 
reactive site which could be a dangling bond, surface void or some defect. The surface gets 
smoothened due to the migration of planar SiHs, as a result the film quality improves. 

B. Formation of bigger clusters 

We have shown in the silane fragmentation process that the product radicals are SiH2, 
SiH2"'", SiHs"*" and SiHs". These radicals are reactive as they have unsaturated bonds and 
can react with silane and with each other forming new species ,0,|2l|. The new 

species formed in the clustering process of the radicals with silane are 

SiH2 + SiH^^ Si2HQ{-2.29eV) (5) 

StH^ + SiH^^ Si2H+{-l.Q8eV) (6) 

SiH+ + SiHi Si2H+ + H2{-0.33eV) (7) 

SiH^ + SiHi Si2H^ + H2{-0.29eV) (8) 
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As we see from the above equations, the interaction of SiH2 with SiH4 resulting in Si2H6 
is highly favored which agrees well with the experimental resultsj^, 0]. 
The species formed by clustering of the radicals with each other are 



SiH2 + SiH2 Si2H,i{-2meV) (9) 

StH2 + SiH^ Si2H+{-3.71eV) (10) 

SiH2 + SiHt Si2H+{-2.90eV) (11) 

SiH2 + SiH^ Si2H^{-2.8QeV) (12) 

SiH^ + SiH^ Si2He{-WmeV) (13) 



The number in the parentheses is the cluster formation energy for the process, which 
is calculated by taking the total energy difference between the product and the reactant 
clusters. 

Except Si2H6 and Si2H6^, other products formed in the above reactions undergo fur- 
ther clustering processes as they have unsaturated bonds. Since Si2H6'^ is saturated with 
hydrogen, it favors fragmentation. Si2H6 is the immediate stable cluster formed by the 
reaction of silane fragmented radicals. It is interesting to see that recent experimental 
observations 

EBB 

|67[ show sufficient amount of Si2H6 in silane decomposed PECVD 

plasma. 



VI. CONCLUSIONS 



Using the ab initio electronic structure methods we have found the ground state struc- 
tures of neutral and charged silane, disilane and their decomposed species. Interestingly, 
in addition to usual Si-H bond, hydrogen forms Si-H-Si and Si-H-H bonds in some of these 
clusters; it forms Si-H-Si bridge bond in Si2H, Si2H2 and Si2H3 clusters and three center 
two electron Si-H-H bond in SiH4"*" and SiHs"*" clusters. We find that charging the clusters 
induces structural modifications which mainly occur through Si-H bond orientations. We 
attribute these structural modifications to electrostatic repulsion between the non-bonded 
electrons and bonded electrons. 

We find that Si-H bond of SimH„ cluster is shortened when positively charged and ex- 
panded when negatively charged. On the other hand Si-Si and Si-H-Si bonds of the cluster 

13 



expand irrespective of the charged state. The bond length variations of Si-H and Si-Si bonds 
in the charged clusters are similar to that of charged SiH and Si2 dimers. 

Our vibrational analysis shows that frequencies of Si-H, Si-Si, Si-H-Si stretch modes show 
significant changes when the clusters are charged. These frequency shifts are consistent with 
their bond length variations. 

We have shown that the low energy pathways of fragmenting neutral and charged silane 
are associated with SiH2, SiH2"'", SiHs^ and SiHs^ products. The fragmentation results indi- 
cate that the production of SiHs can be controlled by electron and H content in the PECVD 
plasma. We speculate that the experimentally observed improvement in a-S:H film qual- 
ity by SiHa deposition is due to the surface migration of planar SiB.^. The investigation on 
the clustering process shows that Si2H6 is the immediate stable cluster formed by the silane 
decomposed radicals which is consistent with the experimental observations^ 
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Table. I. Our calculated ionization potentials ( IP ) of SimH„ clusters using Gaussian98 
package. The experimental values (a) 4^, (b) 41 1 and (c)j4^ are given in the last column. 





Gaussian98 


Experiment 




IP ( eV ) 


IP ( eV ) 


Si 


8.1 


8.2(a) 


SiH 


8.0 


7.9(a) 


SiHs 


9.0 


9(a) 


SiHg 


8.1 


8(a) 


SiH4 


11.0 


11(a), 12.4(b) 


Si2 


7.8 


7.7(c) 


SisH 


8.4 


7.9(c) 


SisHa 


8.0 


7.7(c) 


SisHs 


7.8 


7.8(c) 


Si2H4 


7.9 


7.6(c) 




7.7 


7.7(c) 




9.6 


7.9(c) 
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Table. II The calculated vibrational frequencies v ( cm~^ ) and force constants k ( 
mDyne/A ) of neutral and charged SiH^( m.—1-A) and also charged SiHs clusters using 
B3LYP functional employing 6-31 Ig** basis set. Important modes such as stretch, breathing 
etc., have been identified. 







SiH 


SiH+ 


SiH- 




mode 


V 


k 




k 


u 


k 


1 


Si-H stretch 


2013 


2.49 


2137 


2.81 


1827 


2.0512 






SiH2 


SiH2+ 


SiH2- 


1 


H-Si-H breathing 


1024 


0.64 


907 


0.51 


967 


0.58 


2 


Si-H stretch 


2036 


2.55 


2140 


2.76 


1843 


2.09 


3 


Si-H stretch 


2038 


2.55 


2225 


3.09 


1852 


2.10 






SiHa 


SiH3+ 


SiHg- 


1 


Umbrella 


760 


0.37 


839 


1.11 


894 


0.51 


2 




936 


0.54 


940 


1.11 


972 


0.51 


3 




936 


0.54 


940 


1.05 


972 


0.57 


4 


Si-H stretch 


2192 


2.87 


2274 


3.07 


1885 


2.19 


5 


Si-H stretch 


2233 


3.10 


2356 


3.47 


1885 


2.19 


6 


Si-H stretch 


2233 


3.10 


2357 


3.47 


1899 


2.19 






SiH4 


SiH4+ 


SiH4- 


1 




922 


0.55 


244 


0.04 


712 


0.04 


2 




922 


0.55 


622 


0.26 


796 


0.39 


3 




922 


0.55 


675 


0.28 


814 


0.42 


4 




980 


0.57 


751 


0.34 


992 


0.63 


5 




980 


0.57 


874 


0.48 


1038 


0.64 


6 




2235 


2.97 


1043 


0.65 


1341 


1.11 


7 




2243 


3.12 


2216 


2.96 


1448 


1.25 


8 




2244 


3.12 


2292 


3.29 


2007 


2.47 


9 




2244 


3.12 


3794 


8.55 


2008 


2.46 






SiH5+ 


SiHg- 






1 




63 


0.00 


513 


0.16 
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2 




636 


0.25 


516 


0.16 






3 




656 


0.26 


992 


0.66 






4 




745 


0.37 


1041 


0.71 






5 




869 


0.46 


1041 


0.71 






6 




923 


0.53 


1204 


0.86 






7 




926 


0.53 


1205 


0.86 






8 




1046 


0.65 


1354 


1.09 






9 




2289 


3.12 


1506 


1.38 






10 




2350 


3.45 


1924 


2.29 






11 




2360 


3.48 


1924 


2.29 






12 




3955 


9.29 


1955 


2.27 
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Table. Ill The calculated vibrational frequencies v ( cm~^ ) and force constants k ( 
mDyne/A ) of neutral and charged Si2H„( n=l-4) clusters using B3LYP functional employing 
6-311g** basis set. Si-H-Si ( || ) mode corresponds to the motion of H atom parallel to Si-Si 
bond of the cluster and Si-H-Si ( _L ) corresponds to the motion of H atom perpendicular to 
Si-Si bond. 







Si2 


Si2+ 


Si2- 




mode 


V 


k 


V 


k 


V 


k 


1 


Si-Si stretch 


600 


5.94 


546 


4.92 


586 


5.65 






SiaH 


Si2H+ 


Si2H- 


1 


Si-Si stretcli 


543 


4.02 


407 


2.22 


556 


4.29 


2 


Si-H-Si ( II ) 


1108 


0.76 


890 


0.49 


1104 


0.75 


3 


Si-H-Si ( ± ) 


1565 


1.49 


1384 


1.17 


1479 


1.33 








Si2H2 + 


Si2H2 


1 


Si-Si stretch 


522 


1.93 


395 


1.00 


346 


0.40 


2 




963 


0.59 


993 


0.63 


782 


0.44 


3 


Si-H-Si (II) 


1083 


0.72 


1054 


0.68 


948 


0.55 


4 


Si-H-Si ( II ) 


1175 


0.86 


1171 


0.86 


954 


0.56 


5 


Si-H-Si ( ± ) 


1534 


1.43 


1546 


1.45 


1418 


1.22 


6 


Si-H-Si ( ± ) 


1616 


1.57 


1626 


1.59 


1550 


1.44 






SisHs 


Si2H3+ 


Si2H3- 


1 




358. 


0.07 


198 


0.02 


496 


0.15 


2 


Si-Si stretch 


446. 


1.57 


442 


0.38 


415 


1.53 


3 




647 


0.26 


575 


0.27 


688 


0.29 


4 




739. 


0.35 


637 


0.27 


788 


0.39 


5 




741. 


0.34 


772 


0.38 


802 


0.41 


6 


Si-H-Si ( II ) 


1090 


0.71 


1082 


0.70 


1118 


0.75 


7 


Si-H-Si ( ± ) 


1497. 


1.37 


1606 


1.57 


1375 


1.16 


8 


Si-H stretch 


2082 


2.66 


2192 


2.96 


1888 


2.19 


9 


Si-H stretch 


2089. 


2.68 


2199 


2.97 


1908 


2.24 






Si2H4 


Si2H4+ 


Si2H4 



22 



1 




323 


0.09 


336 


0.08 


198 


0.02 


2 




348 


0.07 


345 


0.07 


360 


0.08 


3 




448 


0.13 


370 


0.08 


458 


0.12 


4 




525 


0.16 


524 


0.17 


676 


0.31 


5 


Si-Si stretch 


563 


0.78 


505 


1.30 


390 


0.95 


6 




616 


0.27 


607 


0.26 


677 


0.30 


7 




919 


0.53 


853 


0.45 


936 


0.54 


8 




956 


0.58 


942 


0.56 


943 


0.55 


9 


Si-H stretch 


2224 


3.00 


2279 


3.14 


1998 


2.43 


10 


Si-H stretch 


2228 


3.01 


2281 


3.14 


2000 


2.47 


11 


Si-H stretch 


2245 


3.13 


2326 


3.38 


2009 


2.46 


12 


Si-H stretch 


2257 


3.17 


2334 


3.40 


2019 


2.52 








Si2H5 + 




1 




128 


0.01 


40 


0.00 


175 


0.02 


2 




390 


0.09 


340 


0.11 


366 


0.70 


3 




406 


0.12 


362 


0.08 


390 


0.09 


4 


Si-Si stretch 


424 


0.41 


397 


0.21 


398 


0.10 


5 




597 


0.25 


608 


0.26 


676 


0.32 


6 




637 


0.28 


676 


0.32 


696 


0.32 


7 




876 


0.49 


817 


0.43 


923 


0.55 


8 




936 


0.55 


910 


0.50 


939 


0.54 


9 




949 


0.56 


921 


0.52 


970 


0.58 


10 




951 


0.55 


976 


0.60 


974 


0.58 


11 


Si-H stretch 


2188 


2.90 


2235 


3.01 


1932 


2.30 


12 


Si-H stretch 


2199 


2.95 


2252 


3.08 


1935 


2.28 


13 


Si-H stretch 


2214 


3.04 


2280 


3.19 


2040 


2.56 


14 


Si-H stretch 


2223 


3.01 


2291 


3.26 


2072 


2.65 


15 


Si-H stretch 


2231 


3.09 


2304 


3.30 


2083 


2.63 








Si2H6+ 





23 



1 




131 


0.01 


105 


0.01 


162 


0.02 


2 




380 


0.09 


208 


0.22 


211 


0.23 


3 




380 


0.09 


265 


0.04 


250 


0.04 


4 


Si-Si stretch 


423 


0.75 


267 


0.04 


285 


0.05 


5 




636 


0.29 


328 


0.07 


498 


0.18 


6 




636 


0.29 


328 


0.07 


562 


0.22 


7 




855 


0.47 


723 


0.34 


855 


0.47 


8 




928 


0.56 


789 


0.40 


884 


0.49 


9 




945 


0.54 


918 


0.52 


932 


0.53 


10 




945 


0.54 


918 


0.52 


938 


0.53 


11 




959 


0.56 


928 


0.53 


962 


0.58 


12 




959 


0.56 


928 


0.53 


966 


0.56 


13 


Si-H stretch 


2209 


2.93 


2256 


3.03 


1765 


1.91 


14 


Si-H stretch 


2217 


2.99 


2258 


3.03 


1826 


2.03 


15 


Si-H stretch 


2218 


3.02 


2314 


3.33 


2058 


2.62 


16 


Si-H stretch 


2219 


3.03 


2315 


3.34 


2061 


2.58 


17 


Si-H stretch 


2227 


3.08 


2321 


3.36 


2069 


2.60 


18 


Si-H stretch 


2229 


3.08 


2322 


3.36 


2076 


2.66 
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FIGURE CAPTIONS 



Fig. 1. The ground state structures of neutral and charged SiH^ ( m = 2-4) clusters. 
Also the structures of SiHs"*" and SiH5~ clusters are shown. All the structures have been 
obtained using Gaussian98 program employing B3LYP/6-311G**. The gray ball represents 
Si atom and small black ball represents H atom. Note that Si and H atoms are connected by 
bonds. In addition to position of atoms, the non-bonded electron density distribution is also 
shown schematically by unconnected half , open and filled circles for neutral, positive and 
negative clusters respectively. The unconnected open and filled circles respectively represent 
the decrease and increase in the non-bonded electron density. 

Fig. 2. Same as in Fig. 1 except for neutral and charged Si2H„ ( n = 1-6 ) clusters. 

Fig. 3. Valance electron density distribution of neutral and charged SiH2 clusters. The 

iso-contour value of 0.07 is chosen to show clearly the density distribution variation between 

the neutral and charged SiH2. The locations of Si and H atoms are also indicated. 

Fig. 4. ( a ) Si-H, ( b ) Si-H-Si and ( c ) Si-Si bond lengths ( A ) of neutral and charged 

Si„H^(n=l or 2, m = 0-6 ) clusters as a function of cluster size. Note that for Si-H-Si bond, 

it is Si-H distance which is plotted.The circles connected by thick hue represents the bond 

lengths of neutral clusters, the square connected by broken line for positively charged and 

triangles connected by doted line for negatively charged clusters. 

Fig. 5. Iso-contour distribution of HOMO in Si2H4, Si2H3 and Si2H3"'". 

Fig. 6. The total energy of SiH3, SiHg and SiHg as a function of Si distance from the center 

of inversion. For each cluster the total energy is given with respect to its ground state energy. 

All these set of calculations have been performed by varying Si distance along symmetry axis 

and freezing H atoms. Since H atoms are frozen in the present set of calculations, for each 

cluster the total energy at the inversion center shown in the figure is slightly overestimated 

in comparison with the fully geometry optimized energy referred in the text. Note that O.OA 

means that Si atom is at the inversion center and the structure of the cluster is planar. 
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